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1.  nmoDDcnoN 


Accurate  predictions  of  tlie  stresses*  natural  Ireqnencies  and  mode  shapes 
of  rotating  blades  for  compressors  are  essential  to  the  performance  of 
aircraft  engines.  The  first  objective  of  this  project  is  to  demonstrate 
that  a  preprocessor  for  the  modeling  of  compressor  blades  with  a  wedge  type 
root  attachment  is  feasible.  A  nser-friendly  program  was  devised  to  be 
used  with  the  minimum  possible  input  that  describes  the  geometry  of 
different  components. 

The  second  objective  was  to  use  this  program  in  conjunction  with  BLADE  to 
perform  steady  stress  and  natural  frequency  analyses  for  a  selected 
example.  The  BLADE  program*  which  is  developed  by  STI*  is  a  special 
purpose  computer  code  for  steam  turbine  blade  stress  and  fatigue  analysis. 
This  program  is  sponsored  by  the  Electric  Power  Research  Institute  (EPRI) . 
A  brief  descriptioc  of  BLADE  is  included  in  the  Appendix. 

The  results  of  these  analyses  were  evaluated  by  their  comparison  to  the 
ANSIS  finite  element  program  results.  These  results  were  validated  using 
experimental  data  provided  by  the  Wright  Research  and  Development  Center  at 
WPAFB. 
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fudib  hjobit 


A  finite  elenent  model  fox  tbe  blede  and  disk  attacbment  is  constructed 
with  siininnsi  input  information  on  the  part  of  the  user.  In  this  section* 
the  user  input  is  described.  An  ezaaiple  which  describes  how  the  finite 
element  model  was  generated  is  given. 

2.1  User  Innnt 

The  user  should  provide  the  program  with  some  information  about  the 
following: 

*  the  vane  sections  geometry 

*  the  root  features 

*  the  disk  geometry 

*  the  blade  materials 

*  forcing 

The  input  file  is  a  free  format  ASCII  file.  It  consists  of  several 
records.  Each  record  starts  with  a  keyword  which  describes  the  input 
section  for  which  the  next  set  of  data  belongs  to.  Tbe  first  record  of  the 
input  file  is  an  exception:  it  contains  the  title  of  the  analysis*  without 
being  preceded  by  a  keyword.  In  the  following*  tbe  input  file  records  axe 
described: 

Record  1  :  TllLB 

Record  2  :  O0?ER  -  Keyword 

IC  -  Material  reference  number 

R1  -  Radios  of  blade  tip  trailing  edge 

R2  -  Radius  of  blade  tip  leading  edge 

Record  3  :  AIRFOIL  -  Keyword  (Figure  2.1.1) 

lA  -  Material  reference  number 

R  -  Radios  of  airfoil  section 

ILB*  ZLB  -  Coordinates  of  leading  edge 

m*  ZIB  -  Coordinates  of  trailing  edge 
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NAP  -  Number  of  defined  airfoil  points  on  the  pressure 

side 

T(l),  Z(l) 

«  • 

•  a 

T(2)f  Z(2).  Pressure  side  coordinates 

KNAP).  Z(NAP) 

NAS  -  Number  of  defined  airfoil  points  on  the 

suction  side 

Kl).  Z(l) 

•  • 

•  • 

«  • 

K2).  Z(2)  Suction  side  coordinates 

KNA8).  Z(NA8) 

Note:  NAP  and  NAS  do  not  have  to  be  equal. 

Record  3  may  be  repeated  for  different  airfoil  sections  up  to  a 
maximum  of  10  sections. 


Record  4  :  ILAXPOBH 
IP 

■D 

■D 

AL 

OB 

cn 

DC 

OFIBT 

RP 

ANB 

ZBID 


-  Kejnrord  (Figure  2.1.2) 

-  Material  reference  number 

-  Radius  to  the  underside  of  platform 

-  Height  of  platform  at  the  downstream  end 

-  Height  of  platform  at  the  upstream  end 

-  Axial  length  of  platform 

-  Distance  from  root  centerline  to  leading  edge  of 
platform 

-  Distance  from  root  centerline  to  trailing  edge  of 
platform 

-  Distance  from  platform  leading  edge  to  the 
coordinate  axis  in  the  axial  direction 

-  Offset  of  root  centerline  from  the  blade's 
global  coordinate  system 

-  Fillet  radius  of  curvature  between  the  airfoil 
and  platform 

-  Slant  angle 

-  The  number  zero  (this  variable  would  be  needed  for 
curved  roots) 
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Figure  I.  .1.1;  Typical  Airfoil  Section  Input 


Record  5  :  BOOT 


miB 

IB 

XI.  n... 

•  • 

•  • 

«  • 

Xt,T8 
H.  82. 

B3.  B4 
VS 

NP 

XU).  Z(l) 

B  % 

•  • 

I(IV),Z(W) 


Keyword  (Figure  2.1.2  and  2.1.3) 
Root  type  nuaiber  (set  to  61) 
Material  reference  number 
Coordinates  of  input  points 


Input  radii 
Stop  width 

Number  of  points  defining  root  end  sections 

Coordinates  of  end  section 


Record  6  :  DISC 
ID 

mv 


Keyword  (Figure  2.1.4) 

Material  reference  number 

Number  of  pair  of  points  describing  the  disk 


Xl  (1).  Z1  (1).  X2  (1).  Z2  (1) 
ZKNFD)  .Zl(lVD)  .X2(HR»)  ,Z2(I1ID) 


Coordinates  of  disk  points 


Record  7 


Record  7 
material 


BATBBIAL  -  Keyword 

II  -  Material  user  reference  number 

IB  -  Reference  number  from  material  library 

n  -  Temperature  at  which  material  is  defined 

may  be  repeated  for  different  material  types  up  to  a  maximum  of  10 
types. 


Record  8  :  SEABE 

FABABBUBS  -  Keyword 

mum 

IS  -  Total  number  of  blades  in  the  stage 

mB 

-  Number  of  blades  per  group 

NBP 

-  Number  of  groups  to  be  analyzed 

m 

-  Number  of  upstream  nozzles 

DBVB 

-  Damping  ratio 

5 


Side 


Figure  2.1.3:  Root  Front  View 


IB 


Rotor  speed 

Blade  of  interest  for  back  substitution 


Record  9  :  FOBCDB 
C 
A1 
A2 

A3 

M 


-  Keyword 

-  Tf  C  =  0.0*  tangential  stream  forces  at  hub  and 
Tip*  respectively:  otherwise*  power  of  the 
stage  and  pressure  drop*  respectively 

-  If  C  =  0.0*  axial  stream  forces  at  hub  and  tip, 
respectively 


Record  10:  EKCnmON  - 
NP 

> 

FIBQ  (1)  - 

mi),  mi) 


Keyword 

Number  of  frequencies  to  be  analyzed  for  dynamic 
stresses 

First  frequency  and  stimulus  ratios  in  the 
Tangential  and  axial  direction 


FBBB  (HP) 
ST  (af).  8A  (mf) 


Last  frequency  and  stimulus  ratios 


Record  11:  BID 


Keyword 


In  the  case  of  natural  frequency  analysis*  records  9  and  10  could  be 
eliminated.  They  can  also  be  eliminated  in  the  case  of  steady  stress 
analysis  under  the  effect  of  centrifugal  forces  only. 

Zs2 _ BmiibIr 

The  blade  to  be  analyzed  is  the  4th  compressor  stage  of  the  FlOO  engine.  A 
typical  airfoil  profile  is  given  in  Figure  2.1.1.  Root  section  information 
is  presented  in  Figures  2.2.1  and  2.2.2.  Figure  2.2.3  shows  the  disk  data 
and  the  input  file  is  given  in  Table  2.2.1. 

The  finite  element  model  was  generated  using  the  preprocessor  written  by 
Sn  for  this  specific  purpose.  The  finite  element  meshes  for  the  blade  as 
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well  88  for  the  different  components  are  given  in  Figures  2.2.4  through 
2.2.10.  Table  2.2.2  gives  the  number  of  eight-noded  isoparametric  finite 
elements  used  in  modeling  each  component. 
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Table  2.2.1,:  Example  Input  File 


AIR  FORCE  COMPRESSOR  GLADE 

COVER 

0 

1 1 .  639  1 1 .  729 


) 

r 


AIRFOIL 

I 

1 I . 292 

-O.  502  -O.  394 

■►O.  510  0.  437 

10 


-0.  500 

-0.  386 

-0.  389 

-0.  300 

-0. 

286 

-0. 220 

-0.  172 

-0.  130 

-0.  072 

-0.  050 

+0.  025 

0  030 

+0. 

145 

0.  130 

+0.  264 

0.  230 

■*■0.  380 

0.  330 

+0.  501 

0.  437 

to 

-0.  495 

-0.  397 

-0.  379 

-0.  330 

-0. 

268 

-0.  260 

-0.  150 

-0.  180 

-0.  041 

-0.  100 

+0.  070 

-0.  010 

+0. 

174 

0.  080 

+0.  290 

0.  190 

+0.  430 

0.  330 

+0  509 

0.  430 

AIRFOIL 


1 

iO.  757 
-0.  500 

-0.  445 

+0.  502 

0.  480 

10 

-0.  499 

-0.  437 

-0.  397 

~0.  340 

-0.  292 

-0.  240 

-0.  198 

-0.  150 

-0.  112 

-0.  070 

0.  000 

0.  034 

+0.  125 

0.  150 

■*■0.  245 

0.  260 

•♦■0.  366 

0.  370 

+0.  490 

0  .  480 

10 

-0.  491 

-0.  447 

-0.  383 

-0.  380 

-0.  268 

-0.  300 

-0.  162 

-0.  220 

“0.  052 

-0.  130 

+0.  050 

~0.  039 

+0.  163 

0.  070 

+0.  276 

0.  190 

♦0.  404 

0.  340 

+0.  5025 

0.  469 

- 

AIRFOIL 

1 

10.  383 

-0.  497 

-0.  479 

+0.  487 

0.  517 

10 

-*0.  496 

-0.  469 

-0.  430 

-0.  398 

-0.  300 

-0.  259 

-0.  205 

-0.  160 

-0.  118 

-0.  070 

-0.  020 

0.  030 

+0.  089 

0.  140 

+0.  210 

0.  260 

■►0  332 

0.  380 

+0.  473 

0.  516 

10 

-0  486 

-0.  480 

-0.  359 

-0.  400 

-0.  240 

-0.  310 

-0.  141 

-0.  230 

— 

-0.  039 

-0  140 

♦0.  064 

-0.  040 

♦0.  184 

0.  090 

+0  292 

0.  220 

■►0.  397 

0.  360 

■»-0,  488 

0.  500 

-lA- 


Table  2.2.1:  Example  Input  File  (Con't) 


AIRFOIL 

I 

9.  848 

-0  492  -0,  524 

■*■0  451  O.  579 
10 


-0 

490 

-0.  512 

-0  430 

-0. 440 

~0.  320 

-0.  320 

-0.  230 

-0.  200 

-0 

130 

-O.  080 

-0.  026 

0,  044 

+0.  090 

0.  180 

+0.  220 

0.  330 

-*■0. 

320 

0.  445 

+0  435 

0.  577 

10 

-0 

480 

-0,  521 

-0.  360 

-0  449 

-0.  244 

-0.  360 

-0.  140 

-0.  270 

-0 

030 

-0.  163 

+0  080 

-0.  041 

+0.  200 

0.  112 

+0.  310 

0.  280 

+0 

389 

0.  420 

-*-0  456 

0.  563 

airfoil 

1 

9  527 

-0.  491  -0.  546 

♦0.  423  0.  612 

10 


-0.  491 

-0.  533 

-0,  399 

-0. 

420 

-0.  301 

-0.  300 

-0.  196 

-0. 

170 

-0  091 

-0.  040 

+0.  004 

0. 

080 

-*■0.  106 

0.  210 

■*■0.  206 

0. 

340 

*0.  325 

0.  500 

+0,  406 

0. 

609 

10 

-0  480 

'•0.  546 

-0.  360 

-0. 

470 

~0.  233 

-0.  370 

-0.  120 

-0. 

260 

-0.  000 

-0.  142 

+0.  107 

-0. 

010 

-*•0.  21 1 

0.  140 

+0  292 

0. 

280 

+0.  370 

0.  440 

+0.  429 

0. 

596 

AIRFOIL 

1 

9.  046 

-0.  484  -0.  568 

+0.  365  0.  656 

10 


-0.  484 

-0. 

554 

-0 

370 

-0.  417 

-0.  291 

-0  320 

-0.  188 

-0  190 

-0  093 

-0. 

060 

+0 

004 

0.  080 

+0.  096 

0.  220 

+0.  183 

0.  360 

+0.  264 

0. 

500 

+0. 

345 

0.  650 

10 

-0.  471 

-0. 

569 

-0. 

330 

-0.  480 

-O.  207 

-0.  380 

-0.  095 

-0  270 

+0  010 

-0. 

150 

*■0 

117 

0.  000 

+0.  208 

0.  160 

+0.  279 

0.  320 

+0  330 

0. 

470 

*■0 

370 

0  640 
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Table  2. 2.1:  Example  Input  File  (Cen't) 


AIRFOIL 

i 

8.  725 

-0.  474  -O,  588 

*0.  291  0.  695 

10 


-0.  473 

-0 

571 

-0  377 

-0.  460 

-0.  265 

-0.  350 

-0.  190 

-O  225 

-0  091 

~0 

OBO 

-0  000 

0  070 

*&.  081 

0.  220 

'*'0.  162 

0,  390 

+0  233 

0 

540 

+0.  276 

0.  686 

10 

-0  460 

-0 

536 

-0.  335 

-0.  510 

-0.  210 

-0.  420 

-0.  105 

-0.  310 

■►0.  011 

-0. 

170 

•*■0.  113 

-0.  010 

+0.  194 

0  160 

+0.  257 

0.  350 

♦O.  291 

0. 

530 

+0.  302 

0.  682 

PLATFORM 

1 

8.372  0.050  0.050  1.402  0.  2764  0  2764  0.634  O  0.172 

21.  0  0 


ROOT 


61  1 

0. 

3925 

0.  2764 

0.  3425 

0.  2764 

0.  3030 

0.  1660 

0.  2625 

0. 

1660 

0. 

1595 

0.  2908 

0.  0806 

0.  2908 

0.  0000 

0.  2287 

-0.  1065 

0. 

1226 

0. 

0300 

0.  0300 

0.  0560 

0.  0560 

0. 

y 

0820 

f 

0. 

OOOO 

1.  1935 

0.  2105 

1.  2701 

0.  2268 

1 .  2760 

0.  2850 

1 

3105 

0. 

3000 

1 .  3252 

0.  3425 

1 .  4020 

O.  3925 

1 .  4020 

OISC 


7. 

1 

5950 

10 

0.  5595 

7.  5950  -0.  5520 

7.  4000 

0.  3574 

7.  4000  -0.  5520 

7. 

3292 

0.  2796 

7.  3292  -0.  2980 

7.  1200 

0.  1515 

7.  1200  -0.  1435 

6. 

9516 

0.  1515 

6.  9516  -0.  1435 

6.  6946 

0.  2351 

6.  6946  -0.  2271 

6. 

5500 

0.  3534 

6.  5500  -0.  3320 

6,  2457 

0.  5530 

6.  2457  -0.  5520 

5. 

8500 

0.  5530 

5.  8500  -0.  5520 

5.  5000 

0.  5530 

5.  5000  -0.  2020 

MATERIAL 

1 

3707 

300 

STAGE 

52  1  0  60 

0.  002  6000  1 

Er-4D 


BLADE 


|ississ:iSHUV.;*A\v4^^ 

fMHBMHHIHIIimtB  \vSK< 

iiSBBMasiiaaaigg 
issitSsiinainHHf'; 


-17- 


Figure  2.2.6:  Finite  Element  Model  of  Root  and  Disk  Section 
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Figure  2.2.7;  Finite  Element  Model  for  Root 
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Airfoil  and  Platform 


Table  2.2.2:  Details  of  Coapresaor  Blade  Model 


Coaq^onent 

Nnaiber  of 

Airfoil 

440 

Platfom 

80 

Root 

726 

Disk 

594 

Total 

1840 

24 
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BLAOB  ANAUrSlS 


The  blade  underwent  two  types  of  analysis,  namely  steady  stress  and  natural 
frequencies  analyses.  Finite  element  programs  ANSYS  end  BLADE  were 
utilized  to  perforc<  the  computations.  Both  programs  produced  identical 
results  which  validate  the  BLADE  program. 

iai _ SteadY.Stress.  Analysis 

The  blade  was  analyzed  statically  under  the  effect  of  centrifugal  forces. 
These  forces  result  from  the  rotation  of  the  rotor.  A  speed  of  6000  rpm 
was  used  in  the  computation.  Steady  stress  variations  are  shown  in  Figures 
3.1.1  through  3.1.3.  1‘aximum  stress  values  are  given  in  Table  3.1.1. 
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Table  3.1.1  :  Maxitrum  Equivalent  Stresses 


Element  U 

Location 

SUcs.^cLl. 

psi 

BLADE 

ANSYS 

1358 

Disk 

10085 

10085 

1359 

Disk 

8604.5 

8604.5 

1546 

Disk 

8460.0 

8460.0 

1361 

Disk 

8411.1 

8411.1 

1362 

Disk 

8375.8 

8375.8 

1360 

Disk 

8315.4 

8315.4 

Maximum  equivalent  stress  in: 

BLADE 

ANSYS 

Airfoil 

4364.9  psi 

4364.9  psi 

Platform 

2530.9  psi 

2530.9  psi 

Root 

7569.0  psi 

7569.0  psi 

29 


hi _ Natural  Freaueiicv  Anelysis 


Natural  Frequency  Analysis  was  peiforoed  on  the  blade  model.  The  equation 
of  motion  for  free  vibratioii  is  of  the  eigenvalue  type, 

[M]{X)  +  [K]{X}  =  0 

where  [N]  and  [K]  are  the  mass  and  stiffness  matrices  and  {X}  is  the 
displacement  vector.  By  solving  this  eigenvalue  problem,  the  eigenvectors 
are  the  mode  shapes  while  the  eigenvalues  yield  the  natural  frequencies. 
Due  to  the  large  number  of  degrees  of  freedom  describing  the  model,  Guyan 
reduction  technique  was  implemented.  At  operation  speed,  centrifugal 
stresses  stiffen  the  blade  and  therefore  the  natural  frequencies  at  speed 
become  higher  than  their  corresponding  values  at  zero  rpm.  The  frequency 
values  at  speed  and  at  zero  rpm  are  given  in  Table  3.2.1.  Mode  shapes  are 
presented  in  Figures  3.2.1  through  3.2.5.  A  Campbell  diagram  which  shows 
the  effect  of  speed  on  the  values  of  the  natural  frequencies  is  given  in 
Figure  3.2.6. 
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Table  3.?..]  :  Natural  Frequencies  for  a  Single  Blade 


- Freonencv  (hz> 


AN  SYS 

BLADE 

at  0  rpm 

at  0  rpm 

at  6000  rpn 

J. 

611.0 

611.0 

657.6 

2 

2173.8 

2173.8 

2212.8 

3 

2610.6 

2610.6 

2625.8 

4 

2854.3 

2854.8 

2866.6 

5 

4939.8 

4939.8 

4969.9 

6 

5466.2 

5466.2 

5502.3 

7 

8565.8 

8565.8 

8597.4 

8 

9313.4 

9313.4 

9347.1 

9 

10247.0 

10247.0 

10257.6 

10 

11118.0 

11118.0 

11141.9 

11 

12865.0 

12865.0 

12902.0 

31 


32 


EiSOd 


-34- 


36 


P0ST4 


Figure  3.2.6:  Campbell  Diagram 


3.3  Comparison  with  Test  Resolts 


Test  results  for  natural  frequencies  test  data  aere  provided  by  AFWAL  in 
the  form  of  Holographic  photos.  The  BI^DE  program  was  utilized  to  plot  the 
displacement  contours  for  the  different  mede  shapes.  Comparison  of  the 
calculated  natural  frequencies  to  the  test  data  is  given  in  Table  3.3.1. 
For  Titanium  material,  it  is  hnown  that  the  elastic  modulus  is  strongly 
dependent  on  the  thermal— mechanical  processing.  For  bar  stoch.  elastic 
modulus  is  about  16X10^  psi.  For  forged  blade  vane,  the  modulus  of 
elasticity  could  become  18.0X10^  psi.  As  a  result,  the  blade  natural 
frequency  can  vary  as  ivch  as  5  percent.  As  a  matter  of  fact,  the  modulus 
for  Titanium  changes  for  one  blade  froD'  one  component  to  the  next  as  shown 
in  Table  3.3.2  obtained  from  tests  conducted  by  Westinghonse.  The  analysis 
was  carried  out  using  a  modulus  of  17.0X10^  psi  as  shown  in  Table  3.3.1  in 
column  (2).  They  compare  favorably  with  the  test  data  as  shown  by  the 
percentage  deviation  in  column  (3).  The  table  also  shows  that  mode  4.  7. 
and  9  were  missing  from  the  experimental  data.  The  Holograms  are  next 
compared  to  the  calculated  displacement  contours  in  Figures  3.3.1  through 
Figures  3.3.8. 


38 


Table  3.3.1:  Comparison  of  Calculated  Frequencies 

and  Test  Data 


Test  Frequency 
hz 

Calculated  Frequency 
bz 

%  Deviation 

602 

611.0 

+1.5 

2138 

2173.8 

+1.67 

2701 

2610.6 

-3.35 

missed 

2854.3 

— 

5200 

4939.8 

-5.0 

5536 

5466.2 

-1.26 

missed 

8565.8 

— 

9480 

9313.4 

-1.76 

missed 

10247.0 

— 

10800 

11118.0 

+2.90 

12620 

12865.0 

+1.94 
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Table  3.3.2: 

Dynamic  (lodnlus 
TI-6AL-4V  Blades 

of  Prototype 

Row 

Nunber 

LocBtion 

Direction 

Average 

Dynamic  Tensile  Modulus  (ksi) 

1 

Airfoil 

long 

17.61 

1 

Root 

long 

17.33 

1 

Root 

trans. 

38.10 

2 

Airfoil 

long 

17.38 

2 

Airfoil 

trans. 

17.73 

2 

Root 

trans. 

17.97 

AO 
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3.3.1:  Comparison  of  Test  Holograms  and  Analytical  Results  for  Mode 


Figure  3.3.2:  Comparison  of  Test  Holograms  and  Analytical  Results  for  Mode 


Comparison  of  Test  Holograms  and  Analytical  Results  for  Mode 


Figure  3.3.4:  Comparison  of  Test  Holograms  and  Analytical  Results  for  Mode 


Figure  3-  '-6:  Comparison  of  Test  Holograms  and  Analytical  Results  for  Mode 


figure  ^.''.7:  Comparison  of  Test  Holograms  and  Analytical  Results  for  Mode  10 


